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In this talk I have raised the question of how the future discovery of leptonic 

\ CP violation can be made robust even at accepting the rather large current ex- 

' perimcntal uncertainties in our knowledges of neutrino propagation in matter. 

\ To make progress toward answering the difficult question, I listed ways to pro- 

■ ceed: (1) Obtain tighter constraints on the MSW theory by testing it by various 

' neutrino experiments. (2) Measure the matter effect in situ, namely within the 

f~i ■ experiment for discovering CP violation itself. (3) Uncover leptonic CP viola- 

tion in a matter effect free environment. I also reported a step made toward the 

^ • above point (2) by taking neutrino factory as a concrete setting; An accurate in 

O I situ measurement of the matter effect looks promising. 
^ '. 

1. Measuring matter density for discovering CP violation? 

Q"^ ■ Leptonic CP violation may be one of the most important elements of our fun- 

O ■ damental understanding o f ma tter. Most probably it is related to baryon number 

^ ■ asymmetry in our universe . Therefore, its exploration is one of the most serious 

1/-^ . targets of future neutrino experiments. 1 argued in the last year in Venice® that lep- 

O I tonic CP violation should exist both in the form of Kobayashi-Maskawa type ® and 

Q I Majorana type ^ phases. It is my pleasure to be here again, by the kind invitation 

j>! I by Milla, to further pursue the line of thought. 1 know that you may say that the 
title of this talk does not sound right for the topics of CP violation. But, in fact, the 



' motivation of my talk addressing this topics does come from the question of how CP 

^ violation (of Kobayashi-Maskawa type) can be uncovered experimentally. 

Let me start by describing my motivation, a personal one. In a recent work 
we have described a concrete proposal for measuring the unknown parameters in the 
lepton mixing and determining the neutrino mass pattern by placing two identical 
megaton-class detectors one in Kamioka and the other in Korea which receive an in- 
tense neutrino beam from J- PARC facility. It may be called as the Tokai-to-Kamioka- 
Korea setting, which is sometimes dubbed as "T2KK" We have demonstrated that 
by using the setting one can resolve the eight-fol d par ameter degeneracy, not only the 
one related to the mass hierarchy and CP phase ^^^^^ but also the 623 degeneracy , 



^ Written version of a talk presented at the "Xll-th International Workshop on Neutrino Tele- 
scope", Venice, Italy, 6-9, March 2007. 



IM^lst^ ^ within reach by the next generation accelerator^ * ' and reactor experiments 

But, one day, I asked myself; "Is this way of uncovering leptonic CP violation 
robust?". Robustness that I concern implies robustness against possible change in 
the basic framework that may arise due to the lack of our knowledge. Of course, full 
treatment of the whole arbitrariness allowed at this moment is just too far a goal, if 
not impossible, to reach. For example, if the three neutrino mixing is too tight to 
accommodate what is happening in nature we stuck, or more appropriately, we face 
with the situation that a gigantic number of completely different scenarios of how CP 
is violated are possible. 

Therefore, I have to be more specific about the setting of the problem. What I 
mean robustness in this talk refers to the one to uncertainties in our underst anding 
of theory describing neutrino propagation in matter, the MSW theory H ^l^^ l'. More 
specifically, I concern uncertainties in the MSW coefficient cmswi which is defined 
by replacing the standard matter potential a = y^GpN^ by cmsw^-- They may arise 
through any effects which renormalize the MSW coefficient. An unimaginative list 
for such effect includes: (1) uncertainty in matter density measurement in the earth, 
(2) presence of non-standard neutrino interactions -. I warn you that the list might 
be short simply because of our ignorance. It is quite possible that there is a real 
candidate for (2) if a new physics at TeV scale is waiting for discovery. 

If the effects outside the standard electroweak theory are absent, looking for the 
deviation of cmsw from unity can serve as a pure test of the MSW theorylB It would 
be desirable if we could distinguish these two aspects, in and outside the electroweak 
theory, in testing the MSW theory. But, in general their effects mix with each other 
and both act as ambiguities in estimating the background to CP violation search. If 
non-standard neutrino interactions have richer flavor structure it might be possible 
to separate between these two effects. 

I believe that it is the very relevant issue not only in T2KK but also in many 
other approaches which seek to uncover CP violation. It is because leptonic CP 
violation is severely contaminated by the matter effect. (For early references on the 
CP phase-matter effect interplay, see e.g., ttSHH /) jf i\^q MSW theory is in error, 
or neutrinos have non-standard flavor changing neutral current interactions, many 
discussions on how to separate genuine CP violation due to the CP violating phase 
from the fake one by the matter effect, most probably, ruin. The potential fragile 
feature of the method for discovering leptonic CP violation is related to the absence 

^ One may argue (as some people in the audience did) that the MSW theory cannot be in error 
because it is based on standard electroweak theory. While it is very likely to be the case, I want to 
emphasize that it is important to confirm our understanding of the theory in the region of coherent 
forward scattering. Recall that despite people believe in Einstein's theory of general relativity, yet, 
they still continue to perform various experimental tests to confirm the correctness of the theory. 
Let me remind you that we never saw a spoon is bent by neutrinos. Or, more precisely speaking we 
saw the spoon may be bent, but only at angle resolution with ^100% error. 



of model-independent (or framework independent) measure for CP violation, lepton 
analogue of "i^'i 27r", an unmistakable clean signature which cannot be masked 
by competing fake effectsE 

You may say that "Oh, the MSW theory is already verified by the solar neutrino 
experiments". Then, I ask "In what sense and to what accuracy?". Let me first 
discuss this question before entering into my problem. 

2. Is MSW theory verified by solar neutrino observation? 

After decades of struggle in solar neutrino observation pioneered by Davis the 
solar neutrino problem is now solved. Among other solar neutrino experiments ^221^ 
SNO finally confirmed that solar neutrinos experience flavor transformation by its in 
situ measurement of CC/NC ratio^^S^ ^j^g phenomenon first discussed by Maki, Naka- 
gawa and Sakata^^H Then, KamLAND pinn ed down the nature of the phenomena as 
due to mass-induced neutrino oscillation . It verified that the phenomenon, which 
was first discovered in the atmospheric neutrino observation by Super-Kamiokande 
exists also with the solar Am^ scale. 

The resultant solution, so called the large mixing angle MSW solution, has a 
unique feature, which is important in our context. The mixing parameters selected out 
live deep inside the adiabatic region. Then, the electron neutrino survival probability 
can be written under the approximation of small ^13 as 



where 6 and 6^ stand for the mixing angle in vacuum and in the solar matter, respec- 
tively. If we further assume that neutrinos are produced in much deeper region than 
resonance density in the sun, we can make further approximation 6m — 7r/2. Then, 
Pee = sin^ 9. That is, the survival probability can be expressed only by the vacuum 
parameter]^ Of course, there is a matter density dependent correction. While it is 
not so small, but it does not appear to be detected by the current experiments. 

Therefore, Pee at high-energy region of neutrinos depends only weakly to the 
matter density in the sun. It is a good news and at the same time a bad news. It is 
a good news because the theoretical prediction of the neutrino rate depends in a 
sensitive manner neither on the absolute matter density in the sun, nor on details of 
the matter density profile in the sun, the celebrated robustness of the prediction by the 

One can argue that T violat ion measurement can provide such clean signature because it is 
not obscured by the matter effect I^Hl While it is in principle true, the experimental setup which 
is required to embody the clean feature of T violation measurement is very demanding. At this 
moment no concrete proposal for such setup is available. 

The fact that Pee can be written only by the vacuum parameter should not be misunderstood 
as absence of the matter effect. Rather, it is the "matter effect dominated" situation in which Pee 
can become much less than 1/2. 



Pee = 7: + 7: COS 26 COS 26. 




LMA MSW solution. On the other hand, it means that observation of neutrinos 
may not be used as a sensitive tool for an accurate test of the MSW mechanism. 
This is the basic reason why the solar neutrino observation can constrain the MSW 
coefficient only up to a factor of 2 uncertainty ^23 _ (Recall that the MSW coefficient 
cmsw is defined by replacing the standard matter potential a = v^G^iVe by cmswo--) 

While it is an independent question, I want to recall that the LMA MSW so- 
lar neutrino solution have not yet been confirmed in a manner independent of the 
standard solar model. Its characteristic predictions, the spectral upturn of the '^B 
neutrinos at low energies and the day-night variation of neutrino fiux of about 2%, 
despite people's great effort, have never been seenB 

To summarize, I have to conclude at this stage that, despite the strong evidence 
for the presence of matter effect in solar neutrino observation, the MSW theory in its 
current status is not established experimentally at the sufficient accuracy to be used 
for reliable estimation of the background matter effect in future CP violation search. 

3. What to be done? 

What I want is simply to obtain robust evidence for leptonic CP violation. By 
"robust" I mean "remain valid even after fully taking account of any experimental 
uncertainties that still exist in the framework I use to define the measurement of CP 
violation". I know that it is not an easy goal to make. But, I believe that discovery 
of CP violation can become a truly experimental statement only when it is done. 

Let me try to describe some considerations to make progress toward the goal. Let 
me restrict myself into the experimental uncertainty in verifying the MSW theory. 
One can proceed along one of the following ways: 

• Prove or obtain much tighter constraints on the MSW theory by testing it by 
various neutrino experiments. 

• Measure the matter effect felt by neutrinos in situ, namely within the experi- 
ment for discovering CP violation. 

• Uncover leptonic CP violation in a matter effect free environment. 

The last possibility sidesteps the problem of uncertainty in experimental verifi- 
cation of the MSW theory in uncovering CP violation; It is certainly a good way 
to proceed. The project which is best suited for this purpose would be the CERN- 
MEMPHYS project El. (The next closest one may be T2K 11^.) Personally, 
however, it is not my choice because of possible drawback of this approach; It is un- 
likely that the same experiment can determine the mass hierarchy (unless one utilizes 

° If one trusts relative normalization between *B and low energy pp plus ""Be neutrino fluxes of 
order ~ 10'' predicted by the standard solar model, then one can show that there is an evidence for 
spectrum upturn by comparing between the SK-SNO and the Ga data. 



the alternative channels, such as high-statistics atmospheric neutrinos or supernova 
neutrinos, etc.). Another relatively matter effect free way of detecting CP violation 
would be the accelerator- reactor combined method ^^HJ 

Unfortunately, I am not able to solve the problem that I raised. The purpose of 
my talk today is merely to pose it in a correct way with expectation that people can 
make progress toward solving the issue. 

4. Neutrino factory measurement of CP violation 

The question I raised might sound too hard to solve. But, the situation is not 
that bad. Let me try to illustrate this point by taking a concrete example. 

The problem of matter effect contamination in uncovering CP violation is most 
serious in neutrino factory ^Si. In a standard setting of placing detector at 3000- 
4000 km from a muon storage ring^^, as is well known, the matter effect dominates 
over the CP phase effect. See, e.g.. Fig. 8 in ^^1^ for illustration of this fact by the 
bi-probability plot. Despite the undesirable feature, people coined into the neutrino 
factory because muon detection is extremely clean, in particular at high energiesE If 
6'i3 is extremely small, sin^ 2^13 <^ 0.01, it is the leading candidate for the machine 
entitled to search for the "diamond" in the frontier. An alternative approach is based 
on "beta beam" EH. 

So far, the problem of ambiguity to CP phase measurement due to the matter 
effect has been addressed in a limited sense, in a form of uncertainty in measured 
matter density by the geophysical method. See e.g., I 9 | 34 | 35 | 36 | 37 | 38 | Here, the 
problem I address is much more broad; In addition to the error of geophysical matter 
density measurement, it includes uncertainty in the theoretical framework in a range 
that is still unconstrained by the current (or the future) experimental measurement. 

5. In situ measurement of the MSW coefficient by neutrino factory 

The fact that the problem is most severe in neutrino factory is, in a sense, a "good 
news" because then one has to necessarily solve the problem. Otherwise, one cannot 
make the goal of unambiguous demonstration of CP violation. In fact, we recently 
made a concrete proposal to solve the problem in situ in measurement in neutrino 
factory . Let me describe our proposal in this section. 

5.1. Which baseline? 

We have started from a general question at which baseline distance the matter 

^ An alternative strategy discussed at the similar stage of understanding how to measure CP 
violation was low-energy superbeam . These two approaches are contrasted in ^ 



density can be measured most accurately|3 To answer the question it is natural to 
consider the response to change in energy because the ratio of the matter effect to 
the vacuum effect, aL/A^i = (^^GpNeL^ / is proportional to neutrino 

energy E. If we measure number of events N[E) and N[E + AE) at two different 
energies E and E + AE, we obtain the double ratio [AN/N) / {AE/ E) where AA^ = 
N{E + AE) — N{E). By taking the double ratio most of the systematic error is likely 
to cancel. The aL dependence of the double ratio (AN/N) / (AE/E) is presented in 
Fig. [H We refer for details of computation. 




Figure 1: aL dependence of the double ratio {AN/N) / {AE/E) defined in the text. Taken from 
EJ 

It is clear from Fig. [1] that the sensitivity to density change is highest at around 
the value of aL where the slope is largest, which occurs at around aL n. It is 
nothing but the one called as the "magic baseline" in the literature^ Interestingly, 
the magic baseline appears in our context the bast distance for measuring the earth 
matter density, which is not surprising because it was known to be the characteristic 
length of the matter effect 

5.2. How can the matter density he measured accurately? 



^ When we talk about the measurement of matter density by neutrinos what we mean is, of 
course, the electron number density in the earth. Since the electron fraction is very close to 1/2 
in the earth it can be relate d to the matter density. 

^ It has been proposed EESl lY^aX a second detector at the baseline can serve as a powerful 
degeneracy solver because of its special property of independence of CP phase 5 . The distance 
has been known in the theory of neutrino propagation in rnatter as "refraction length" ^^'^ . For a 
recent discussion on the meaning of the magic baseline, see SSI. 




Figure 2: The energy distribution of event number (per GeV) is plotted with three values of the 
matter density, 4.2 g/cm"^ (shown in blue dotted curve), 4.3 g/cm^ (black solid), and 4.4 g/cve? 
(red dash-dotted). The left and the right two panels in Fig. [2] are for the cases of the normal 
and the inverted mass hiera rch ies, respectively. The mixing parameters are taken as (5 = and 
sin^ 26li3 = 0.01. Taken from SSI 

The next question we must address is how the matter density can be measured 
accurately by neutrino factory. The answer to this question can be found in Fig. Ej 
in which the energy distribution of event number (per GeV) is plotted with three 
values of the matter density, p = 4.2g/cm (shown in blue dotted curve), 4.3g/cm 
(black solid curve), and 4.4g/cm^ (red dash-dotted curve). It is indicated in the 
figure that the matter density dependence of the oscillation probability changes at 
a critical neutrino energy Ec of ~ 10 GeV and ~ 20 GeV in the neutrino and the 
anti-neutrino channels, respectively. At E > Ec, higher the matter density, smaller 
(larger) the oscillation probability in neutrino (anti-neutrino) channel. At E < E^ 
the behavior is reversed; higher the matter density, larger (smaller) the oscillation 
probability in neutrino (anti-neutrino) channel. These are the case of normal mass 
hierarchy. In the case of inverted mass hierarchy, the above described behavior is 
completely reversed as seen in Fig. [2J It was shown in that these behavior can be 
simply understood by using the approximate analytic formula derived in for the 
Ue appearance probability. 

From Fig. [2] the appropriate analysis principle is obvious; two energy bin analysis 
with neutrino and anti-neutrino running combined. Because of the opposite response 
to the matter density change in the neutrino and the anti-neutrino channels the 



analysis will lead to a compact allowed region in sin^ 2^13 — p plane, which is confirmed 
by the actual analysis. See Fig. 2 of 

People may ask the question; Now, we are taking the different method of analysis 
from the one considered in the previous subsection, the energy scan. Then, is the 
distance comparable to the magic baseline still the bast place for accurate measure- 
ment of the earth matter density? The answer to this question is provided in Fig. [31 
In both cases with and without varying 6 the best sensitivity to Ap/p is achieved at 
baseline L = 7500 — 9000 km. It confirms our expectation, but note that there is 
nothing sacred in the magic baseline. 
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Figure 3: Plotted are the L dependence of Ap/p, the figure taken from xhe red solid line is 
for the case in which S is fixed to be 100 degree. The blue dotted line is for the case in which d is 
varied with weight of gaussian distribution centered at the above value with width of 20 degree 
Sni. Xhe normal hierarchy is assumed and ^13 is taken as sin^ 2^13 = 0.01. 

5.3. Analysis results 

Let me jump to the analysis results, leaving the details of the procedure to the 
description inl^. In our analysis the sensitivity to 5 possessed by the near detector 
placed at 3000-4000 km is modeled by adding the a gaussian error for 6 with width 
of 20 degree. In Fig. H] presented are the fractional errors Sp/p of the matter density 
determination as a function of sin^ 26'i3, with three curves corresponding to Icr, 2a, 
and 3a CL defined with 1 DOF by marginalizing 61^. They correspond, roughly 
speaking, the best cases in each mass hierarchy. 

We notice that determination of the matter density p in very long baseline neutrino 
factory represented in Fig.|l]is extremely good; The uncertainty 6p/ p is only about 1% 
level even at 3a CL at the largest value of ^13, sin^ 26'i3 = 0.1. The uncertainty remains 
small, about 1% at sin^ 2^13 = 0.01 at la CL for both the normal and the inverted 
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Figure 4: The fractional errors in the matter density determination Sp/p at 1, 2, and 3 a CL 
defined with 1 DOF by marginahzing 613 are plotted as a function of sin^ 2^13 by the red sohd, the 
green dash-dotted, and the blue dotted lines, respectively. The upper panel is for the normal mass 
hierarchy with S — and the lower panel for the inverted mass hierarchy with S = 47r/3. 

mass hierarchies. At sin^ 2^13 = 0.001, the uncertainty increases to about 2% (2.5%) 
at 1(7 CL for the case of the normal (inverted) mass hierarchy. At sin2 2^i3 = 0.0001, 
however, Sp/p becomes worse to about 3% (4%) at la CL for the respective mass 
hierarchies, which however is still within a tolerable level for CP analysis. 

5.4- Unexpected 6 dependence at the magic baseline 

Unfortunately, it is not the end of the story. Look at Fig. O in which a curious 
dependence of 6p/p on the CP phase 6 is reported; 6p/p blows up to a rather large 
value at some particular region of 6. The strong 6 dependence looks like "against the 
definition" of the magic baseline. 

The reason for such curious behavior is, however, understandable. Notice first 
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Figure 5: Presented are the fractional errors Sp/p at 1 a CL with 1 DOF as a function of 6 
for five different values of ^13, sin^ 2^13 = O.f (red solid line), sin'^ 29i3 — O.Of (green dashed line), 
sin^ 2^13 = 0.003 (blue short-dashed line), sin^ 26*13 = 0.001 (magenta dotted hue), sin^ 26*13 = 0.0001 
(light-blue dash-dotted line). The upper and the lower panels are for the normal and the inverted 
mass hierarchies, respectively. 

that the "disease" occurs only at small values of ^13, sin^ 2^13 < 0.003. At such 
small ^13, every term in the appearance oscillation probability is small, and more 
specifically, the dominant atmospheric term and the solar-atmospheric interference 
term are comparable. Therefore, the response to matter density change of a particular 
term can be cancelled by that of the other term, producing insensitivity to the matter 
density change. It is not surprising that this phenomenon occurs at some particular 
values of 6, which depend upon the mass hierarchy. It is shown in that the 
relationship between the value of 6 at which the "disease" occurs in the normal and 
the inverted mass hierarchies (which differs by ~ vr) can be understood in this way. 

A more sophisticated analysis of matter density measurement in neutrino factory 
is carried out in by explicitly combining yields at near (4000 km) and the far (7500 



; normal hierarchy - 




± 



inverted hierarchy 
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km) detectors. The problem of 6 dependent loss of the sensitivity to 6p/p is certainly 
cured to some extent as can be seen by comparing Fig. [6] to Fig. [5l in agreement with 
the argument given in^^. At the same time, however, we have to conclude that the 
problem is not completely resolved by the improved treatment. 

L=7500km + L=4000km, sin^ 26*13 =0.001 




1 2 3 4 5 6 
True value of 6cp 



Figure 6: Presented are the fractional errors Sp/p at la (yellow region), 2a (green region), and 3a 
CL (blue region) with 1 DOF as a function of S for sin^ 2^13 = 0.001. The normal mass hierarchy 
is assumed, and therefore, the la line should be compared to the magenta dotted line of the upper 
panel of Fig. [5] The results is obtained by using the same analysis procedure as in SSI and the figure 
is by courtesy of Walter Winter. 

I want to emphasize, however, that apart from the 6 dependent "disease" , neutrino 
factory measurement at around the magic baseline is able to resolve in situ the issue 
of uncertainty not only of the earth matter density but also of the MSW coefficient. 
Finally, I have two remarks: 

• I note that there is a valid criticism to our viewpoint I just described above, the 
one we took in^®. That is, we propose to measure the matter density along the 
neutrino trajectory for baseline of 7500 km. But, the sensitivity to CP violation 
in neutrino factory is mainly possessed by the near detector at 3000-4000 km 
from the muon source. Since the neutrinos that are "sensitive to CP violation" 
pass through a different part of the earth, the matter density determined for 
7500 km baseline does not solve the ambiguity issue that arises in search for CP 
violation. 

Now, I certainly agree with this argument itself, in principle. But, I also want 
to point out that the following; The two trajectories for 7500 and 3000-4000 



km baselines differ in the ratio of path lengths in the lower mantle to the upper 
mantle regions. It is believed that the difference between the upper and the 
lower mantle regions arises because of the change in the form of matter, from 
silicate to perovskite that occurs because of the higher pressure in the deeper 
mantle region. If the states of the matters are better understood by ongoing 
study of matter under high pressure (see for example one will be able 
to relate the matter density in the lower mantle region to that in the upper 
mantle region, opening the possibility of relating the matter density measured 
by detector at baseline of 7500 to that along the trajectory of 3000-4000 km. 

• It is important to clarify the connection between the averaged matter density 
along the trajectory in the earth and the "matter density" measured by neu- 
trinos, p = NgYe (assumi ng w e know Ye). It appears that the latter is larger 
than the former by ~ 5%^^. The theoretical understanding of this point is, 
however, not very transparent at this point. 



6. Is T2KK advantageous? 

I have motivated to testing the MSW theory to make future discovery of leptonic 
CP violation a robust experimental evidence, not just an intriguing hint. While it is 
quite a generic problem for every project to search for CP violation, my concern is 
in particular on T2KK ^ as is natural for one of the proponents. The question is 
whether T2KK is the better setting for this purpose. 

I argue here that T2KK is indeed advantageous. In the T2KK setting, the sen- 
sitivity to 5 mainly comes from the Kamioka detector, while the Korean detector is 
indispensable for the mass hierarchy resolution. Because of the relatively short base- 
line (L = 295 km) of the Kamioka detector, the uncertainty of the matter density will 
not produce any serious ambiguities to the CP sensitivity. Moreover, the assumed 
value of the matter density that comes from geophysical estimation together with the 
uncertainty of the MSW coefficient as a whole can be cross checked by an in situ 
measurement of the matter density in T2KK, as we examined for neutrino factory. 
In this task and in detecting CP phase effect, the comparison between the Kamioka- 
Korea two identical detectors will play a decisive role; the (anticipated) power of the 
two-detector method'^. Therefore, it will provide another robust way of identifying 
leptonic CP violation due to the leptonic Kobayashi-Maskawa phase. 

I must admit that my above argument is merely an argument. It does not make 
much sense unless it is backed up by the real quantitative analysis. I hope that I can 



come back to this issue in the next year in Venice. 



7. Conclusion 

In my talk, 1 have raised the question of how the future discovery of leptonic CP 
violation can be made robust even accepting the rather large current experimental 
uncertainties in the matter density and the MSW coefficient. To make progress toward 
the difficult goal I have suggested the several ways to proceed. (1) Obtain tighter 
constraints on the MSW theory by testing it by various neutrino experiments. (2) 
Measure the MSW coefficient in situ, namely within the experiment for discovering CP 
violation itself. (3) Uncover leptonic CP violation in a matter effect free environment. 
I also reported a step made toward the above point (2) by taking neutrino factory as 
a concrete setting. I hope that people warmly accept the legitimate question, take 
it seriously, and can make progress toward the goal of robust demonstration of CP 
violation.® 

Finally, I have to give a cautionary remark; If my discussion gave you the impres- 
sion that the uncertainty in the theory of neutrino propagation in matter is the only 
potentially important obstacle to clean discovery of CP violation in future neutrino 
experiments, it is certainly misleading, and I have to apologize for it. Experimentally, 
more urgent issue would be to control the systematic errors related to neutrino flux 
and cross sections. Fortunately, great amount of efforts are dedicated to improve the 
situation and people are making progress ^ 48.49,50j _ 
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